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A.. EXPERIMENTAL AM) ANALYTICAL 
INVESTIGATION OP AX I SYMMETRIC D1PEUSERS 

L. A. Povlnelii 

lewis Rruirch Centtr 
Cleveland, Ohio 44135 


A finite difference computer program for tur- 
bulent compreaalble flow waa need to eatabllah the 
perf orma-te of aeveral dlffuaer ahapea for experi- 
mental taating. The dlffuaera were dealgned to 
have a linear change in Mach number, a linear change 
In preaaure or a curvature fitted by a quadratic 
equation. Teatlng waa performed with M • 0.1 to 
0.9 with and without boundary layer bleed. Above 
M ■ 0.6, data were obtained with a normal ahock up- 
etream uf the dlffuaer entrance. Peak atatlc prea- 
aure recovery occurred with a dlffuaer inlet M- 
0.75. The quadratic dlffuaer yielded the hlgheat 
total preaaure recovery. 

Introduction 

The efficient perfo-mance of a auperaonlc In- 
let eyetem la dependent on the proper dealgn of the 
aubaonic dlffuaer aa well aa the auperaonlc inlet 
portion. Although a considerable amount of effort 
has been directed toward the auperaonlc Inlet, a 
relatively smaller effort hue been devoted toward 
studying aubaonic dlffuaera for advanced supersonic 
aircraft. Relatively few studies have been con- 
cerned with axlsynmetrlc dlffuaer design. O-J) One 
of those studies has discussed the possibility of 
Improving supersonic Inlet performance by the use 
of short subsonic diffusers. (*) Furthermore, little 
is known regarding the effects of boundary bleed, 
inlet flow distortions, free-atream turbulence and 
vortex generators on axleymmet rlc diffuser perfor- 
mance. 

This study la directed towards arriving at de- 
sign criteria for axlaymuetrlc dlffuaera for ad- 
vanced supersonic aircraft. Specifically, the 
pressure recovery and minimum distortion were meas- 
ured for several diffuser geometries and the results 
compared with the predictions from an vlscid- 
lnvlacld computer program. The effects of boundary 
layer bleed and vortex generators were also meas- 
ured. Tests were conducted over a Mach number range 
of 0.1 to 0.9. Above Mach 0.6, a normal shock waa 
present upstream of the diffuser Inlet. The Rey- 
nolds number fur the tests was from 0.3 X 10* to 
1 x 106. 

Analysis 

Computational Method 

The finite-difference procedure of Anderson^) 
for turbulent, swirling, compressible flow In axl- 
syimsetric ducts was used to compute diffuser per- 
formance. 

The equations of motion are solved In a coor- 
dinate system made up of the streamlines and poten- 
tial lines of the lnvlscld solution. Boundary lay- 
er type approximations are made In that coordinate 
system; the viscous effects being treated as per- 
turbations on the flow. The analysis solves for the 
entire flow across the duct at each streamwlse ata- 
* 

Aerospace Research Engineer, Associate Fellow AIAA. 


tlon. This strong interaction solution eliminates 
the matching problems resulting from the lnvlscld 
flow merging with the boundary layer. 

Diffuser Geometries 

Performance calculations have been carried out 
for a series of axisymmetrlc diffusers and are pre- 
sented In Ref. 6. Four different types of diffusers 
were chosen for Investigation; namely a linear Mach 
number change, a linear area change, a linear pres- 
sure variation and a quadratic area change. Only 
two of the four types analysed in Ref. 6 were select- 
ed for this study. The first type selected was the 
linear Mach number change, or dM/dz diffuser, be- 
cause it Is representative of current designs. The 
oM/dz diffuser is characterized by gradual curvature 
at the entrance and low Initial rates by diffusion. 
Maximum curvature and diffusion rate occur well down- 
stream of the entrance, and the skin friction under- 
goes a relatively gradual decrease. The second type 
of diffuser selected for study was the quadratic ar- 
ea variation, that Is, the wall contour waa estab- 
lished by requiring the area distribution along the 
diffuser to be specified by a quadratic equation. 

This type of diffuser, which differs considerably 
from the dM/dz, Is characterized by rapid diffusion 
and large w.ll curvature at the entrance and la sim- 
ilar to that studied by Stratford. The skin 
friction dropa very rapidly to a low but finite val- 
ue, resulting In minimum energy dissipation. The 
proximity to zero skin friction coefficient makes 
the diffus. r sensitive to flow dlsturbsnces or dis- 
tortions and may cause separation. However, the 
ability to successfully operate a diffuser of this 
type would lead to a maximum amount of diffusion In 
a minimum "ength. High recoveries could, therefore, 
be accomplished In a shorter length than other con- 
ventionally shaped diffusers. 

In this report the per* rmance of the dM/dz 
and tne Stratford (quadratic area variation) dif- 
fusers will be presented. 

A pparatus and Procedure 
General Approach 

Testing was carried out using the apparatus 
shown in Fig. 1. Entering air was accelerated by 
the converging-diverging nozzle to supersonic veloc- 
ity. The exhaust pressure was regulated so that a 
normal shock occurred In the divergent portion of 
the nozzle The air then entered Into the diffuser; 
where both the Inlet and exhaust properties were 
measured in order to obtain performance data for the 
dllfuaer only. 

Test Geometries 

The diffuser geometry consisted of a straight 
cowl (12” dlam) and Interchangeable center bodies. 
Four centerbody shapes were tested and are shown in 
Fig. 2. The coordinates for the two diffusers are 
given In table I. The ratios of the diffuser length 
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to the annular height (at the diffueer entrance) 
were 6 and 12. The entrance annular height wee 
conatant (II) in.) (or all the geoswtrlea. The 
range of length to height ratloa epane the range of 
eubeonic axleyasetric dlffuaera teated on euperson- 
ic lnleta at the NASA- Lewie and Ames Research ten- 
ters. The dllfuaer exit-to-entrence area ratio was 
2.0 (or all geometries. 

The removeable centerbody was preceded by a 
ehort (2.3" long) boundary layer bleed section and 
the convergent-divergent nozsle (see Tig. 1). The 
bleed region had 600 hnlea drilled norsial to the 
eurface. The holes had an 0.067” diameter and were 
arranged in a staggered 2 row pattern. The amount 
of bleed was variable from sero to 3 percent of the 
inlet flow. The slope of the bleed section was 
3 degrees. 

Test Conditions 

Air was supplied at a maximum total pressure 
of 60 palg (33 lbs/sec) to the diffuser. The exit 
flow was removed by a 26” Hg vacuum (maximum) ex- 
hauster. The Reynolds number (based on annular 
height) was varied from 3 x 10* to 1 x 10®. The 
range of Mach numbers at the diffuse! entrance was 
varied from 0.1 to 0.9 (see Measurements section). 
Below a diffuser entrance Mach number of approxi- 
mately 0.6 the flow was subsonic throughout the ap- 
paratus. Above 0.6, a normal shock (or shock train) 
was present in the divergent section of the noirle 
or In the bleed region. The diffueer entrance Mach 
number was Increased by reducing the back pressure. 
The reduction in back pressure also caused the shock 
location to move further downstream. At the high- 
est Mach number used in this testing, the shock oc- 
curred in the bleed region. It Is noted that the 
diffuser pressure recovery Is defined in terms of 
the diffueer entrance and exit total pressures and 
does not take into account the pressure loss across 
the shock wave. Therefore, higher recovery does 
not necesearlly mean higher pressure recovery at the 
compressor; particularly above Mach 0.623. 

Measurements 

A eight tube rake at the diffuser entrance 
was used to measure seven total pressures equally 
spaced across the inlet annulus and the static 
pressure at the centerline of the annulus. Wall 
static pressures were measured every 1 In. along 
the cowl surface, from the nozrle to the diffuser 
exit. Circumferential static pressure meuaurmcnis 
were also made on the cowl at 2, 6, and 10 in. down- 
stream of the diffuser entrance. Wall static meas- 
urements were made on the centerbody as -'.own in 
Fig. 2. Dynamic pressure measurements were made at 
two locations on the centerbody eurface. Circum- 
ferential static pressure measurements were made at 
2 or 3 axial stations (sec Fig. 2) on the diffuser 
centerbody eurface. An exit rake was located either 
at the 9 in. or 18 in. location. The rake was com- 
posed of 16 total pressure tubes distributed across 
the annulus and a static pressure tube at the cen- 
terline of the annulus. Wall static pressures were 
measured on the centerbody and on the cowl at the 
same axial stations corresponding to the inlet and 
exit rakes. A linear interpolation was made with 
the two wall statics and the mid-stream static in 
order to obtain static pressure- at the total pres- 
sure tube locations on the rake. These data were 
then used to determine local Mach numbers across the 
annulus. The local values were then area-weighted 


to yield area-weighted Mach numbers at the dlffuaer 
entrance and exit planes. The data presented in 
this report are mainly plotted as a function of the 
area-weighted-dlff user-entrance Mach number or sim- 
ply the diffuser entrance Mach number. The data ob- 
tained at the Inlet and exit rakee were alec used to 
compute area-weighted pressure coef f iclente, static 
preesure rise and total pressure recovery. 

Vortex Generators 

The diffuser geometries we'e tested with and 
without vortex generators. For '.he short diffuser, 
one row of generators was located approximately 1" 
downstream of the diffuser Inlet. With the long 
diffusers the first row was at 0.8 In. for the dM/dt 
diffuser and at l 1 ; in. for the quadratic geometry. 

The second row was located 7.6 In. downstream of the 
diffuser Inlet. The vortex generators had a chord- 
to semi span rttio of 2. The generator cross sec- 
tion was equivi lent to half of a NACA Series 0012 
airfoil. The leading edge waa rounded to have a 
0.003" to 0.006” R leadlig edge. The generator 
chord length was 1/2 in. and the ratio of vortex 
generator spaclng-to-chord length was nominally 1.7. 
The generators were oriented at 16 degreee angli of 
attack and were in counter-rotating pattern. 

Results 

Static Pressure Distribution 

A typical static pressure distribution through 
a section of the nonle, bleed section and las long 
dM/dz diffuser is shown in Fig. 3. The distribution 
starts just downstream of the con.ergent part of the 
centerbody and extends throughout the divergent noz- 
zle section, the bleed section and the diffuse' as 
indicated on Fig. 3. The diffuser inlet station 
corresponds to an X/L value of 0 and extends to e 
value of 1. Figure 3 shows the distribution obtained 
with a diffuser entrance Mach number of 0.867 aa de- 
termined from the diffuser entrance preesure rake 
(see Measurement section). The pressure distribution 
shows that the flow is accelerated to supersonic ve- 
locities in the divergent section. A normal shock 
(or sh- ck train) In thr bleed region reduces the 
flow to subsonic velocity before entering the dlf- 
fuset (X/L • 0). The static pressure rise then oc- 
curn In a continuous fashion through the diffuser. 

li l et Pressure Profiles 

The total pressure profile measured by the rake 
at the diffuser entrance is shown in Fig. 6(a). The 
inlet profile is shown for three values of the Mach 
number at the diffuser entrance. Distortion of the 
p.oflle Is seen to Increase with Mach number ax the 
.hock system strengthens. This distortion is be- 
lieved to be caused by the presence of a normal 
shock or shocks in the bleed region. Measurement of 
the pressure profile upstream of the nozrle land up- 
stream of the normal shock retlon) showed uniform 
profiles at all the test conditions. The effect of 
using a small amount of boundary layer bleed (1/2 
and 22 of the entrance flow at M 0.8 and M 0.C re- 
spectively) on the centerbody (contoured surface of 
the diffuser) Is shown in Fig. 6(b). There is little 
effect of the bleed flow on the pressure profile at 
M 0.6. However, at M 0.7 the entire cross-stream 
profile is afficted. The boundary layer on the cen- 
terbody e .ace le greatly reduced and the mld- 
strea . profile is more uniform. At M 0.83, the 
amount of bleed available was ”ery small, and the ef- 
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fact on the profile la minimal. Tha raaulta pre- 
aantad In the succeeding aactlona were, therefore, 
obtained with uniform Inlet profllee up to M O.A - 
0.7. At M0./» 0.8, the prof Ilea were very dls- 
torted In the abaencr of boundary layer bleed and at 
M 0. 8 - 0.9 the Inlet flow waa dlatorted with or 
without boundary layer bleed. 


Total Preaaure Recove ry 


The total preaaure recovery for the long 
Stratford dlffuaer (L/AR ■ 12) la ahown In Fig. 5(a), 
Excellent recovery waa obtained from Mach 0.15 to 
0.70 (98 to 99 percent). Above Mach 0.7 the recov- 
ery decteaaea very rapidly. Consular analyaia of 
thla dlffuaer Indicated that the recovery would re- 
main high over the entire Mach range trated, If the 
Inlet profile waa uniform. The entrance flow to 
the dlffuaer, however, becomes very dlatorted at 
the higher Mach nuabera. The drop off In recovery 
la believed due to the aenaltlvlty of the dlffuaer 
to Inlet flow dletortlona. The Stratford dlffuaer 
haa a vary rapid rate of dlffualon Immediately u. 
the entrance. Thla rapid dlffualon coupled with the 
onaet of a normal shock (or ahock train) makea the 
dlffuaer very aenaltlve to flow diatortlona, there- 
by increaalng the auaceptlblllty to flow aeparatlon. 
The effect of boundary layer bleed on recovery la 
ahown in Fig. 5(b). An Improvement In recovery la 
obaerved over the range from M 0.75 to M 0.85. Aa 
aeen In Fig. 4, the bleed flow cauaed the Inlet pro- 
file to become more uniform at M 0. 7 whereas, the 
limited bleed at M 0.85 did not cause much change. 
Hence, the recovery of the dlffuaer, which la very 
aenaltlve to inlet diatortlona, Increased In a cor- 
responding fashion. Installation of a row of vortex 
generators at the upstream location did not improve 
the recovery. The lack of Improvement la due to 
the fact that the generators could not be located 
sufficiently far upstream for the vortices to erzr- 
gize the boundary layer before aeparatlon, whi h 
probably occurred at the sharp corner. 

The pressure recovery for the long dM/dz dif- 
fuser la shown In Fig. 6(a). The same general trend 
la observed In the recovery behavior as with the 
Stratford diffuser. However, the values of the re- 
covery are lower up to a Mach number of 0.75. A 
crossover occurs at that point and the dM/dz dif- 
fuser exhibits better performance than the 
Stratford. Hence, the sensitivity of the dM/dz dif- 
fuser appears to be much less than that of (he 
Stratford diffuser. The absence of high rates of 
dlffualon at the entrance and the gentle Initial 
curvature of the dM/dz geometry make It relatively 
Insensitive to the inlet flow distortions. Indeed, 
the use of bleed did not have l significant Influ- 
ence on recovery. Computer analysis showed that 
the dM/dz diffuser would have higher skin frlitlon 
and hence lower recovery than the Stratford diffus- 
er over the entire Mach range tested. This '.ulavior 
was verified by the results above p to Mach 0.75. 
The analysis for the dM/dz geometry also Indicated 
flow separation at a pol..c near the dl! fuser exit, 
where Che centrrbody has u rapid final curvature. 

One, therefore, vu Id expect vortex generators to 
be effective, ii Vested sufficiently far ups'ronm. 
The experimental results with n row of vortex gen- 
erators in the upstream location Improved the re- 
covery 1/2 to 3/4 percent: (see Fig. 6(b). 

An additional feature waa obaerved in the re- 
covery curves. At a Mach number of approximately 
0.625, a normal shock occurs upstream of the dif- 


fuser Inlet. A alight rise In recovery la observed 
at that condition Thla Increase may result from a 
reduced akin friction coefficient In the region of 
the normal shock. Ms thews at al ,M< have analysed 
the normal shock boundary. layer Interaction data of 
Seddon'" 1 using Coles universal wake function. 

The results (see fig. 1 of ref. 8) show that the In- 
tegrated akin friction coefficient over the region 
where the ahock affeeta the boundary layer la less 
than th. skin friction coefficient with an undisturb- 
ed boundary layer. The reduction In akin friction 
coefficient leads to lower shear atrass at the wall 
and higher total preaaure recovery. 

Pressure Kecoverv - Short Diffusers 

Computer analyaia indicated that flow separa- 
tion would occur close to the entrance of the short 
Stratford dlffuaer (L/AR - 6). Hence, the use of 
boundary layer bleed might he expected to prevent 
aeparatlon. The use of vortex generators, on the 
other hand, would not be expected to be helpful due 
to an Insufficient path length to enable the vor- 
tices to energize the low momentum boundary layer. 

The experimental results in Fig. 7 bear out the ex- 
pectations. The recovery la generally low (Fig. 7 
(a), no bleed, no generators) and la verv sensitive 
to Inlet flow distortions. Boundary layer bleed 
Improved the performance appreciably (Fig. 7(b)) by 
allowing the flow to remain attached at the dlffuaer 
entrance. However, the dlffuaer la still quite sen- 
sitive to diatortlona at the higher Mach numbers. 

The effect of shock-boundary layer interaction on 
akin friction la readily apparent at M 0.625 In 
- g. 7(a). 

The results for the short dM/dz diffuser (L/ 

AH - 6) Is ahown In Fig. 8(a). The pressure recov- 
ery la generally 1 to 2 percent lower than that ob- 
tained with the longer geometry (L/AR * 12). Use of 
the computer analysis Indicated that separation 
would occur near the exit of the dlffuaer (X/L 
0.9). Hence, locating vortex generators near the 
dlffuaer entrance should Introduce sufficient mix- 
ing to prevent separation. The «xperlm«ii:al results 
show that the generators do enhance the performance; 
(see Fig. 8(b). The recovery with generators la 
within 1/4 to 1/2 percent of those obtained with the 
larger dM/dz diffuser (L/AR - 12) over the entire 
Mach number range. The short dM/dz diffuser alao Is 
Insensitive to distortions In the entrance flow. 
Therefore, a reduction In diffuser length by a fac- 
tor of 2 does not appear to degrade the diffuser's 
performance. The rise in pressure recovery, when 
the normal shock occurs In the throat region (M 
0.625), is seen In Fig. 8(a). As mentioned previ- 
ously, this rise la attributed to skin friction re- 
duction. From the results of Matthews 18 ), It ap- 
pears as If skin friction la re 'uced when the bound- 
ary layer la thickened. Use of bleed, therefore, 
should decrease boundary layer thickness, yielding 
higher friction losses and poorer recovery. The ex- 
perimental results wl’h bleed show that, Indeed, the 
recovery is 1/2 to 1 percent lower than the no bleed 
results. It Is further noted that bleed Is generally 
Ineffective with the dM/dz diffuser because the 
gradual curvature and lack of a high dlffualon rate 
allow for the development of a healthy boundary lay- 
er. 


As Indicated earlier (Apparatus and Procedure), 
hlg.ter diffuser recovery does not necessarily mean 
higher pressure at the compressor face. 
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Distortion 


Nomenclature 


The long diffusers yielded distortion vslues 
(l.s., Ifj/Pj) thst were 7 percent st M 0.7, ss 
shown in 71g. 9(e) and (b). The date for the short 
diffusers Is given In Fig. 10(a) for the quadratic 
end 10(b) for the dM/dt. At M 0.7, the distortion 
Is 15 percent for the quadratic and 12 percent for 
the dM/dt diffuser. 

Coapsrlson with Computations! Analysis 

Figure 11(a) shows a comparison of the exper- 
imental results for the dM/dt diffuser and the com- 
puter analysis. The experimentally measured total 
pressure profiles at the diffuser entrance were 
used as Input for the coaiputat Ions. The computer 
results show a rise In recovery st an entrance Mach 
number of approximately 0.62, corresponding to the 
experimental results. The Increase In pressure re- 
covery is believed to be associated with the In- 
creased blockage (thlcxer boundary) layer), leading 
to a reJuctlon In skin friction coefficient. Fig- 
ure 11(b) compares the results for me long quad- 
ratic diffuser. Again, the experimental pressure 
profiles at the diffuser entrance were used as 
starting conditions for the numerical calculations. 
T he experimental results and the computations show 
reasonably good agreement over the Mach number 
range . 
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static pressure 
total pressure 
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Conclusions 

The results from this study have shown thst 
dlf fuser-length-to-helght ratios of 12 tend to be 
very conservative In design. 

The highest total pressure r overy, up to an 
entrance Mach number of 0.75. -a obtained with a 
Stratford diffuser (L/oa - id). At a Mach number of 
0.7, the recovery was 96 percent and the total pres- 
sure distortion at the diffuser exit was 7 percent. 
The long dM/dx diffuser (L/AR • 12) had a slightly 
lower recovery; being 97 percent at Mach 0.7 with s 
distortion of 7 percent. The dM/dx diffuser, how- 
ever, was less sensitive than the Stratford diffus- 
er to the distortions present in the entrance flow 
at the higher Mach numbers (0.75 to 0.9). The In- 
creased sensitivity of the Stratford diffuser Is 
caused by a combination of the rapid divergence at 
the entrance and the effect of the upstream normal 
shock on the behavior of the boundary layer. 

A reduction of 50 percent In the length of dM/ 
dz dltf- er (L/AR - 6) did not cause any significant 
reduction In the performance. At Mach 0.7, the 
pressure recovery was 97 percent; the distortion was 
12 percent and the diffurer was Insensitive to en- 
trance flow distortions. Therefore, for s slight 
Increase in distortion level, one c.m obtain com- 
parable recovery with half the length. The short 
Stratford diCfuser yielded leas favorable results. 
The recovery was the lowest (95S percen' M 0.7) 
and the distortion was the highest (15 pe.,-ent). 

The experimentally obtained pressure profiles 
st the diffuser entrance were used as Input condi- 
tions for a vlscld/lnvlscld diffuser program. Cood 
agreement was obtained between the numerical compu- 
tations and the experimental data for the long dif- 
fusers. Use of the analysis for the short diffus- 
ers indicated where flow sefsratlon would occur; 
thereby providing information for the location of 
vortex generators. 
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Figure 1. * Experimental diffuser apparatus. 
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Figure 6. - Total pressure recovery for long dMfdz diffusers. 
UAR • 12. no bleed. 
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Figure 8. - Total pressure recovery for short dMIdz dif- 
fusers, L/AR • 6, no bleed. 
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Figure 9. - Total pressure distortion at diffuser exit station 
for long diffusers. UAR ■ 12. 
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(b) d/Wdz DIFFUSER. WITH VERTEX GENERATORS, NO 
BLEED. 

Fio'jrc lU. - Total pressure distortion at diffuser exit sta- 
tion for shori diffusers, L/AR ■ 6. 



<a> dM/dz DIFFUSER. NO GlNERATORS. NO BLEED. 
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(b) STRATFORD DIFFUSER. NO VORTEX GENERATORS. 
no Bleed. 
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Figure 11. - Comparison of experimental data and numer- 
ical computation, L/AR ■ 12. 
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